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The MSW solution to the solar neutrino problem in the presence of 
random solar matter density perturbations 
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We present the evolution equation describing MSW conversion, derived in the framework of the Schrodinger 
approach, in the presence of matter density fluctuations. Then we analyse the effect of such fluctuations in the 
MSW scenario as a solution to the solar neutrino problem. It is shown that the non-adiabatic MSW parameter 
region is rather stable (especially in 5m 2 ) for matter density noise at the few percent level. We also discuss the 
possibility to probe solar matter density fluctuations at the future Borexino experiment. 
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1. The present deficit of solar neutrinos seems 
to disfavour any " astrophysical solutions" || 
whereas it points to neutrino oscillations. In par- 
ticular the resonant conversion due to neutrino 
interactions with constituents of the solar mate- 
rial (the Mikheyev-Smirnov-Wolfenstein (MSW) 
effect) H is the most elegant and viable expla- 
nation for the existing solar neutrino data. It 
provides an extremely good data fit in the small 
mixing region with Sm 2 ~ 10 _5 eV 2 and sin 2 29 ~ 

io- 3 -nr 2 

Here we investigate the stability of the MSW 
solution with respect to the possible presence of 
random perturbations in the solar matter density. 

The existence of matter density perturbations 
at the level of 1% or so cannot be excluded ei- 
ther by the Standard Solar Model (SSM), which 
is based on hydrostatic evolution equations, or by 
the present helioseismology observations |ll[] . 

Let us remind that in Ref. the effect of peri- 
odic matter density perturbations added to an av- 
erage density po, i.e. p(r) = po[l + ^sin(7r)] upon 
resonant neutrino conversion was investigated. In 
that case parametric resonance in the neutrino 
conversion can occur when the fixed frequency (7) 
of the perturbation is close to the neutrino oscil- 
lation eigen-frequency and for rather large ampli- 
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tude (h ~ 0.1 — 0.2). There are also a number of 
papers which address similar effects by different 
approaches [5 10 1. 

In the present discussion we consider "white 
noise" matter dentity perturbations Sp (as in Ref. 
[jTcj] ) . Namely we assume that the random field 
Sp has a <5-correlated Gaussian distribution: 



{5 P {r x )8p{T2)) = 2p 2 (e)L S(r 1 - r 2 ) 



(1) 



where £ = Sp/p and the correlation length Lq 
obeys the following relation: 



free 



< L < A n 



(2) 



In (|2|) the lower bound is dictated by the hydro- 
dynamical approximation used later on, Zf ree = 
(cm) -1 being the mean free path of the particles 
in the solar matter On the other hand, the 
upper bound expresses the fact that the scale of 
fluctuations have to be much smaller than the 
characteristic v matter oscillation length A m , as 
indeed the 5-correlated distribution in eq. ([[]) re- 
quires. 

2. According to the standard Schrodinger 
equation approach, we derive now the most gen- 
eral neutrino evolution equation in random mat- 
ter density. The evolution of a system of two 



1 For Coulomb interactions, the cross-section a is deter- 
mined by the classical radius of electron ro e = e 2 /m e c 2 ~ 
2 X 10~ 13 cm, resulting in 1{ ICC ~ 10cm for a solar mean 
density no ~ 10 24 cm -3 . 
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neutrinos v e and v x (x = fi or r)^| in the solar 
matter is governed by 



v. 



H e H ex 
H-ex 7l x 



(3) 



dt \ V, 

where the entries of the Hamiltonian matrix are 

H e = 2[A ex (t)+A ex (t)}, H x = 0, 
5m° 



AE 



sin26>, 



I 8 771^ 

Ae X (t) = -[Ve X {t)-—COs2e] 

A ex (t) = -V ex (t)-£. 



(4) 



Here 9 is the neutrino mixing angle in vacuum, 
8m 2 the mass squared difference, E the neutrino 
energy and the matter potential for the v e — > v x 
conversion reads 



V ex (t) = ^Lp(t)(Y e ), 



(5) 



where m p is the nucleon mass and Y e is the elec- 
tron number per nucleon. 

The above system can be rewritten in terms of 
the following equations: 

j(t) = H e (t)R(t)-2H ex (t)(P(t)- 1/2) 

R(t) = -H e (t)I(t) 

P(t) = 2H ex (t)I(t), (6) 

where P — \ v e \ is the v e survival probability, R = 
Re(^ e ) and 7 = Im(is x V e ), with the corresponding 
initial conditions P(to) = 1> 7{to) — 0, R{to) = 0. 

The Eqs. (|^) have to be averaged (see Q for 
more details) over the random density distribu- 
tion, taking into account that for the random 
component we have: 

(i^ +1 ) = 0, (ies(i)ies(ii)) = 2n5(t - h), (7) 
where the quantity n is given by: 

«(t) = (A 2 ex (t))L Q = \v e l(t)(e)L . 



(8) 



2 Here for simplicity we consider only the case of solar v e 
conversion into active state or v T \ however the discus- 
sion can be extended also to the case of conversion into a 
sterile state 



In terms of the averaged quantities defined as 

<p(t)) = r(u), (R(t)) = n(u), (/(*)> = z(u), 

we can write the variant of the set (|6[) as: 

±{t) = 2[A ex (t)K(lA) - «(U)J(U) - H 1s (P(u) - oo/G)] 
K(t) = -2[A ex (t)l(U) + k(u)TZ(u)\ 
V{t) = 2H ex X(U). (9) 

This system of equations^] explicitly exhibits the 
noise parameter k. It is now possible to outline 
the main effects due to the presence of the ran- 
dom field Sp upon the resonant neutrino conver- 
sion. The MSW resonance condition remains un- 
altered, i.e. A ex (t) = V ex (t) -5m 2 cos 28/ 2E = 0, 
due to the random nature of the matter perturba- 
tions. Due to the condition Lq <C A m , the noise 
parameter k (cfr. Eq.(||)) is always smaller than 
A ex (t) except at the resonance region. As a conse- 
quence, (see Eqs. (||)) the perturbation can show 
its maximal effect just at the resonance provided 
that the corresponding noise length 1/n obeys the 
following adiabaticity condition at the resonance 
layer Ar 



Ar(/t) 



> I. 



(10) 



This condition is analogous to the standard 
MSW adiabaticity condition at resonance a r = 
Ar/(A m ) res > 1 H). For definiteness, we have 
taken Lq = 0.1 x A m . The two adiabaticity pa- 
rameters are related as 



e 



tan 2 29 



6m 2 sin 2 26R 
AttE cos 29 



(11) 



Therefore due to the restriction (||) and for the 
range of parameters we are considering, £ ~ 10 -2 , 
tan 2 29 > 10~ 3 — 10~ 2 , we have a r < a r . 

As a result of (|TT|), in the adiabatic regime 
a r > 1, the effect of the noise is enhanced to the 
extent that the mixing angle is small. Further- 
more, the MSW non-adiabaticity a r < 1 always 
brings to a r < 1. As a result in our discus- 
sion the fluctuations are expected to be ineffec- 
tive in the non-adiabatic MSW regime. Finally, 
it can be shown that the matter noise weakens 



3 These equations are equivalent to those obtained in 
Ref. [p.011 in terms of the variables x = 21Z, y = —TL and 
oo. 
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the MSW suppression in the resonance layer, ex- 
hibiting somehow the role of a friction. 

3. In view of the qualitative features just out- 
lined, we discuss the implications of noisy solar 
matter density in the MSW scenario for the solar 
neutrino problem. We have solved numerically 
the coupled differential Eqs. in (0) for the v e sur- 
vival probability, using as reference SSM the most 
recent Bahcall-Pinsonneault model (BP) ||. 

The x 2 analysis has been performed taking 
the latest average d expe rimental data of chlo- 
gallium 



rine |1 
experiments 
(74±8)SNl 
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14 and Kamiokande JL5J 
(2.55 ± 0.25)SNU, R e ^ = 



i2££ = (0.44±0.06)ig£ 
Fig. 1 
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The results of the fitting in the 5m 2 , sin 2 29 pa- 
rameter space, are shown in Fig. 1, where the 
90% confidence level (C.L.) areas are drawn for 
different values of £. One can observe that the 
small-mixing region is almost stable, with a slight 
shift down of 5m 2 values and a slight shift of 
sin 2 29 towards larger values. 

The large mixing area is also pretty stable, ex- 
hibiting the tendency to shift to smaller 5m 2 and 
sin 2 29. The smaller Sm 2 values compensate for 
the weakening of the MSW suppression due to 



4 For gallium result we have taken the weighted average of 
(77 ± 8 ± 5)SNU|Q and SAGE 



GALLEX datum 



Ga 



(69±11±6)SNU 



the presence of matter noise, so that a larger por- 
tion of the neutrino energy spectrum can be con- 
verted. The £ = 8% case, considered for the sake 
of demonstration, clearly shows that the small 
mixing region is much more stable than the large 
mixing one even for such large value of the noise. 
Moreover the strong selective 7 Be neutrino sup- 
pression, which is the nice feature of the MSW ef- 
fect, is somehow degraded by the presence of mat- 
ter noise. Consequently the longstanding conflict 
between chlorine and Kamiokande data is exac- 
erbated and the data fit gets worse. Indeed, the 
presence of the matter density noise makes the 
data fit a little poorer: Xmin — 0- 1 f° r £ = 0j 
it becomes Xmin = 0.8 for £ = 4% and even 
X 2 mm = 2 for £ =8%. 

In conclusion we have shown that the MSW 
solution exists for any realistic levels of matter 
density noise (£ < 4%). Moreover the MSW so- 
lution is essentially stable in mass (4- 10- 6 cV 2 < 
5m 2 < 10~ 5 eV 2 at 90% CL), whereas the mixing 
appears more sensitive to the level of fluctuations. 

4. Let us also stress the fact that the solar neu- 
trino experiments could be viable tools for provid- 
ing information on the matter fluctations in the 
solar center. In particular, the future Borexino 
experiment jl6), aiming to detect the 7 Be neu- 
trino flux could be sensitive to the presence of so- 
lar matter fluctuations, as the 7 Be neutrinos are 
those mostly affected by the presence of matter 
noise. 

In the relevant MSW paramter region for the 
noiseless case, the Borexino signal cannot be defi- 
nitely predicted (see Fig. 2a) . Within the present 
allowed C.L. regions (dotted line) the expected 

rate, Z Be = R^/RbZ™ ( solid lines )' is in the 
range 0.2 0.7. 

On the other hand, when the matter density 
noise is switched on, e.g. £ = 4% (see Fig. 
2b), the minimal allowed value for Zs e becomes 
higher, Z Be > 0.4. Hence, if the MSW mecha- 
nism is responsable for the solar neutrino deficit 
and Borexino experiment detects a low signal, say 
%Be ^ 0.3 (with good accuracy) this will imply 
that a 4% level of matter fluctuations in the cen- 
tral region of the sun is unlikely . 

Once more, the solar neutrino detection turns 
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Fig. 2a 




sin 20 



Fig. 2b 
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out to be an important approach for studying the 
solar physics. 

We thank Z. Berezhiani, P. Krastev, S. 
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